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Abstract. Quasi-free photoproduction of 7/ mesons off nucleons bound in the deuteron has been measured 
with the combined Crystal Barrel - TAPS detector. The experiment was done at a tagged photon beam 
of the ELSA electron accelerator in Bonn for incident photon energies from the production threshold up 
to 2.5 GeV. The ?/-mesons have been detected in coincidence with recoil protons and recoil neutrons. The 
quasi-free proton data are in good agreement with the results for free protons, indicating that nuclear 
effects have no significant impact. The coincidence with recoil neutrons provides the first data for the 
771 — > nrj reaction. In addition, also first estimates for coherent ^'-production off the deuteron have been 
obtained. In agreement with model predictions, the total cross section for this channel is found to be very 
small, at most at the level of a few nb. The data are compared to model calculations taking into account 
contributions from nucleon resonances and t-channel exchanges. 

PACS. I3.60.Le Meson production - 14.20.Gk Baryon resonances with S=0 - 25.20.Lj Photoproduction 
reactions 



1 Introduction 

The complex structure of the nucleon is still one of the 
greatest challenges for the understanding of the strong in- 
teraction in the low energy, non-perturbative regime. One 
expects that, like in nuclear structure physics, the main 
properties of the interaction are reflected in the excita- 



Correspondence to: B. Krusche, Klingelbergstrasse 82, CH- 
4056 Basel, Switzerland, e-mail: Bernd.Krusche@unibas.ch 



tion spectrum of the nucleon, but so far the correspon- 
dence between model predictions and experimentally ob- 
served states is quite unsatisfactory. All constituent quark 
models predict more states than have been observed. This 
problem of 'missing resonances' becomes more severe the 
higher the excitation energy. However, the experimental 
data base is dominated by elastic scattering of charged 
pions off the nucleon, which profits from large hadronic 
cross sections, but is biased against states that couple only 
weakly to Nn. The combination of continuous wave elec- 
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tron accelerators with sophisticated An detection systems 
now allows the study of photon induced reactions with at 
least comparable precision as hadron induced reactions. 
Therefore, photoproduction of mesons has developed into 
a prime tool for the investigation of the nucleon excitation 
scheme [3(2]. 

Photoproduction of light mesons like pions at high in- 
cident photon energies involves many partial waves, so 
that the interpretation of the data requires sophisticated 
partial wave analyses. Such programs are under way and 
will largely profit from the combination of polarized pho- 
ton beams with polarized targets giving access to sin- 
gle and double polarization observables. However, alter- 
natively due to the suppression of higher partial waves, 
the photoproduction of heavier mesons close to their pro- 
duction thresholds may give access to resonances which 
contribute only weakly to other channels. Photoproduc- 
tion of 77-mesons, which is completely dominated in the 
threshold region by the Sn(1535) resonance, is the best 
studied example for this approach [H][Hl[71Rf^lfTTTlP^lP^l[Tr?l 
[T4l[T5l[T6l[T7l[T8l[T^[2f)l[2T] . Since the mass of the rf (m v > 
958 MeV [3 ) is much higher than the 77-mass (m v ~ 548 
MeV [3]), resonances contributing to rj' threshold produc- 
tion may have masses around 2 GeV. Of course, lower 
lying resonances may also contribute due to their large 
widths. Because of their iso-scalar nature, both, r\ and rf 
offer the additional selectivity that only N* resonances 
can couple to Nrj,Nr]'. Excited Z\-states can emit these 
mesons only when decaying to other A's, in particular the 
4(1232), and thus contribute to the ?/7r-channel but not 
to single r/ -production (again such processes have been 
recently intensively studied for reproduction in the tjtt - 
channel [T4l|2^|2^|2^|231|2^r2T] ) . Therefore, rf threshold 
production is expected to have a large sensitivity to N* 
resonances at excitation energies, where the missing res- 
onance problem is most severe. This is illustrated in Fig. 
[TJ where the experimentally observed nuclear excitation 
scheme is compared to model predictions. 

Until recently, rj' photoproduction was not much ex- 
plored, not even for the proton. In an early attempt, Mukho- 
padhyay and coworkers [28 analyzed bubble chamber data 
with an effective Lagrangian model and concluded that 
the dominant contribution comes from the excitation of a 
Di3(2080) resonance. Analyses of a more recent measure- 
ment with the SAPHIR detector 29,30 claimed contri- 
butions from different resonances (Sn, Pn) and strong t- 
channel contributions. However, these data are not in good 
agreement with three later measurements with the CLAS 
detector at Jlab |3TUT3] and the Crystal Barrel/TAPS 
setup at ELSA [TB]. These second generation experiments, 
which profit from much better counting statistics and bet- 
ter control of systematic effects due to the use of highly ef- 
ficient detector systems with large and uniform solid angle 
coverage, clearly supersede the previous data. Nakayama 
and Haberzettl [35] presented an analysis of the earlier 
CLAS-data in the framework of an effective Lagrangian 
model. They found possible contributions from Sn, Pn, 
P13, and D13 resonances. However, also these results are 
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Fig. 1. Comparison of experimentally established nucleon res- 
onances [3] (left and right) to predictions in the framework of 
non-relativistic quark models [4| (center). The shaded (color: 
green) area indicates the range accessible for rf photoproduc- 
tion in the present experiment. 



far from being unique since the available cross section data 
do not sufficiently constrain them. 

Additional information may be obtained by exploring 
the iso-spin degree of freedom. The electromagnetic ex- 
citations of N* resonances are iso-spin dependent. Reso- 
nances which are only weakly excited for the proton may 
give stronger signals for the neutron and vice versa. Inter- 
ference patterns between different resonances and between 
resonances and background contributions may change due 
to sign changes of the electromagnetic couplings. Again, 
■q photoproduction with the Sn(1535) - Di 3 (1520) inter- 
ference 33 , 34 , 35 and the prominent excitation structure 
above the Sn range, which is only seen for the neutron 
35,36,37 , is a very instructive example. 

Due to the technical problems involved in the measure- 
ment of small production cross sections off bound nucle- 
ons, photoproduction of 7/-mesons off the neutron had not 
been studied up to now. Here, we report the first results 
for quasi-free ?7'-photoproduction off protons and neutrons 
bound in the deuteron. The paper is organized in the fol- 
lowing way. Sections [2] and [3] summarize the experimental 
setup and the data analysis. The bases of the models the 
data are compared to are discussed in section |U The re- 
sults are summarized in section [5] First, the data in co- 
incidence with recoil protons are compared to free proton 
data as a check of the quasi-free production hypothesis. 
Both, quasi-free proton and neutron data are then com- 
pared to model fits. In addition to the quasi-free data, a 
first estimate for the cross section of the coherent process 
7c? — > drj 1 was extracted. Final conclusions are given in 
section [BJ 
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Fig. 2. Overview of the experimental setup at the ELSA accelerator. 



2 Experimental setup 

The experiment was done at the electron stretcher ac- 
celerator ELSA in Bonn [38,39 . For the measurements 
discussed here, electron beam energies of 2.6 GeV and 
3.2 GeV have been used. Real photons were produced 
with the bremsstrahlung technique. Their energies were 
tagged via the momentum analysis of the scattered elec- 
trons by a magnetic spectrometer, which is schematically 
shown in Fig. [3] (see [15] for details). For this experiment 
only the part of the focal plane covered by the scintillating 
fiber detector but not the part covered by the proportional 
wire chamber was used. This limited the maximum tagged 
photon energies to 80 % of the electron beam energy. The 
different beam settings are summarized in Tab. [TJ 



Table 1. Summary of beam times. E e ~: electron beam en- 
ergy, E lt : maximum energy of tagged photons, E po r. energy 
of maximum linear photon beam polarization, <2> : energy inte- 
grated electron flux. Total life time: beam time multiplied by 
acquisition life time. 
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Fig. 3. Setup of the tagging spectrometer. 



Due to these settings and the typical 1/-E 7 intensity 
distribution of the photon flux, the average time inte- 
grated flux at photon energies above 2 GeV was roughly 
lower by a factor of two than the flux at lower energies. 
The largest part of the 3.2 GeV beam time was done with 
a copper radiator foil of 0.3 % radiation lengths thickness, 
producing unpolarized bremsstrahlung. For a small part of 
this beam time and for the running with the 2.6 GeV elec- 
tron energy a diamond radiator was used to produce a lin- 
early polarized photon beam via coherent bremsstrahlung 
(see [30] for details about linear polarization) for the ex- 
traction of photon beam asymmetries. However, the sta- 
tistical quality of this observable was marginal for rjf pro- 
duction, since for most of the beam time the setting of the 
polarization peak was optimized for ^-production at lower 
incident photon energies (see Tab. [TJ). Therefore photon 
asymmetries have not been analyzed for the ?/-channel. 
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The target consisted of a vertically mounted cryostat 
attached to a tube entering the Crystal Barrel detector 
from the upstream side. The target cell itself was a cap- 
ton cylinder (0.125 mm foil thickness) with a diameter of 
3 cm and a length of 5.3 cm, filled with liquid deuterium 
(surface thickness 0.26 nuclei/barn). The reaction prod- 
ucts emerging from the target were detected with electro- 
magnetic calorimeters covering almost the full solid angle; 
the Crystal Barrel (CB) detector (1290 Csl crystals cov- 
ered the full azimuthal angle for polar angles between 30° 
and 168°) gj and the TAPS detector (528 BaF 2 crystals 
mounted as hexagonal forward wall covered polar angles 
down to 4.5°) [421143] . Plastic scintillator detectors in front 
of the TAPS modules and a scintillating fiber detector [4~4~] 
inside the Barrel covering the same polar-angle range were 
used for charged particle identification. A schematic view 
of the full arrangement is shown in Fig. [5J more details can 
be found in [35] , where apart from the target an identical 
setup was used. The time-of-flight walls were mounted but 
not used for this experiment. 




Fig. 4. Front view of the TAPS forward wall. Left hand side: 
logical segmentation for the LED-low trigger, right hand side: 
logical segmentation for the LED-high trigger (see text). 



The first-level hardware trigger for the experiment was 
exclusively based on signals from the TAPS forward wall 
detector. The reason is that the CB was read out by photo- 
diodes without timing information. Measurements of reac- 
tions off the free proton can use signals from recoil protons 
traversing the Inner detector. However, in order to have an 
identical trigger setting for quasi-free production off the 
proton and off the neutron this option was not used. The 
modules of the TAPS detector were equipped with two in- 
dependent leading edge discriminators, combined in two 
different ways into logical groups (see Fig. @|. For most 
of them (ring 12-5 from outer edge to center) the lower 
threshold was set to «55 MeV (LED-low) . It was set to 80 
MeV, 135 MeV, 270 MeV for rings 4, 3, 2, respectively. The 
inner-most ring was not allowed to trigger. The LED-high 
thresholds were set to 70 MeV for rings 9-7, rising from 
105 MeV (ring 6) to 180 MeV (ring 2). Again, the inner- 
most ring was not allowed to trigger and the three outer 
rings (block G) had no LED-high. The first level trigger 
included two components: one or more LED-low discrim- 
inators from at least two logical sections above threshold 
or at least one LED high discriminator above threshold. 
In the second case, a second-level trigger from the FAst 



Cluster Encoder (FACE) of the Crystal Barrel, indicating 
at least two separated hits in the Barrel, was required in 
addition. All first level triggers thus required detection of 
at least one photon in TAPS. Such a trigger is only ef- 
ficient for reactions with a high multiplicity of photons 
like the rf — > ir°ir°ri — > 67 or the 77 — > 3ir° — > 67 decays. 
But even then the trigger efficiency for mesons at back- 
ward angles is not large. In principle, also events where 
a recoil nucleon is detected in TAPS might activate the 
hardware trigger. This would, however, lead to uncontrol- 
lable trigger efficiencies since the LED thresholds are only 
calibrated for photons (recoil nucleons have different sig- 
nal shapes in BaF2 scintillators) and the energy deposited 
by neutrons is more or less random. Therefore only events 
where photons alone (identified by non-firing veto detec- 
tors and time-of-flight versus energy) fulfilled the first level 
trigger conditions were accepted in the analysis. 

3 Data Analysis 

3.1 Particle and reaction identification and extraction 
of cross sections 

Photoproduction of ?/-mesons was identified via the rj — s> 
?77r 7r — > 67 decay chain, which has a branching ratio of 
8%. Cross sections were extracted for four different reac- 
tion types. The two most important ones are quasi-free 
production off the proton 7c? — > (n)prj' , which requires co- 
incident detection of an rj and a recoil proton and quasi- 
free production off the neutron 7c? — > n(p)r]' via detection 
of the rj together with a recoil neutron. For the control 
of systematic uncertainties (see below) also the inclusive 
reaction jd — > {np)r)' with no condition for recoil nucleons 
was analyzed. Finally, also an estimate of the fully inclu- 
sive reaction 7<i — > Xrj was obtained, where also final 
states like tj'tt contribute. 

In the first step of the analysis photon and neutron 
candidates (called 'neutral hits') were separated from pro- 
ton candidates (called 'charged hits'). This was done in the 
CB with the help of the scintillating fiber detector and in 
TAPS with the charged particle veto detectors. In TAPS 
a hit was assigned to 'charged' if the veto of any cluster 
module or the veto of any neighbor of the central module 
of the cluster had responded (even if the neighbor module 
itself had no signal above threshold). The latter condi- 
tion applies to charged particles which traverse the edge 
of a veto but deposit their energy in the neighbor module 
(due to large impact angles). All other hits were assigned 
to 'neutral'. In the Barrel, a hit was assigned to 'charged', 
if at least two layers of the Inner-detector had recorded 
a hit within ±10° of azimuthal angle. It was assigned to 
'neutral' if no layer had fired within this azimuthal angle. 
Hits with one responding layer of the Inner-detector were 
discarded. In the TAPS forward wall, correct identification 
of protons, neutrons, and photons can be additionally con- 
trolled with a time-of-flight versus energy analysis, while 
in the CB no direct separation of photons and neutrons is 
possible (more details on particle identification are given 
in 05] (CB) and gS] (TAPS)). 
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this condition, the 'best' combination was chosen via a x 2 
minimization: 



0.8 0.9 1 1.1 

m(jtVTi)[GeV] 



-0.5 0.5 
Am[GeV] 



-0.5 0.5 
Am[GeV] 



Fig. 5. Top to bottom: reactions ^d — > (np)r)' , 'yd — > p(n)rj' , 
*yd — > n(p)ri' . Left hand side: 67 invariant mass spectra, (Blue) 
dashed curves: background ht, (red) solid histograms: sum of 
background and jy'-peaks, in addition: simulated line-shapes 
(pink). Center: missing mass spectra for cut on rf invariant 
mass peaks, right hand side: (black) points: missing mass spec- 
tra for events in invariant mass peak after background subtrac- 
tion, solid (red) curves: simulated line-shapes. All spectra for 
incident photon energies from threshold to 2 GeV (integrated 
over all beam times and the full polar angle range). 



In the next step, events with at least six 'neutral' hits 
as candidates for the f/-decay photons were selected and 
assigned to four partly overlapping classes corresponding 
to the above defined reaction types: six 'neutral' and one 
'charged' hit for the {n)prf final state, seven 'neutral' for 
the n{p)rj' final state, six or seven neutral or six neutral 
and one charged for (np)r)', and at least six 'neutral' with- 
out any further condition for Xrf . 

The identification of the rfN final states was then 
based on a combined invariant and missing mass analysis. 
The invariant mass analysis identified the rf , the missing 
mass analysis excluded events where further mesons have 
been produced but have escaped detection (except for the 
Xrf final state where such events were included). 

The invariant mass of all possible disjunct photon pairs 
was calculated. Only events having at least one combina- 
tion of six 'neutral' hits to two photon pairs with invariant 
masses between 110 and 160 MeV (pions) and one pair be- 
tween 500 and 600 MeV {rf) were kept. In cases where the 
photons could be combined in more than one way to fulfill 



m fc (7 7 )-M fe ) 2 



3 

V- 

£J (Z\m fe (77)) 2 



(1) 



where for each disjunct combination of the six photons 
into three pairs the invariant masses are ordered so that 
w 1,(77) — ™2(77) < ^3(77). The Z\mfc(77) are their 
uncertainties (computed event- by- event from the detector 
resolution for energies and angles) and Mk is the 7r°-mass 
for k = 1,2 and the 77-mass for k = 3. The above case 
applies to events with exactly six 'neutral' hits, where in 
total 15 different combinations are possible (events with 
recoil proton or without detected recoil nucleon) . For can- 
didates for the quasi-free reaction off the neutron (seven 
'neutral' hits) one must in addition loop over the unpaired 
hit, since in CB photons and neutrons cannot be distin- 
guished. This corresponds at maximum to 105 combina- 
tions, giving rise to larger combinatorial background. In 
this case, the hit which was not assigned to a pion- or 
77-decay photon is accepted as neutron candidate. In the 
case of the Xrf final state even higher multiplicities may 
occur. 

As in |45j the nominal masses of the mesons were then 
used as a constraint to improve the experimental resolu- 
tion by re-calculating the measured photon energies from 



/ 1 
E l2 = Ei t 2- 



(2) 



where E\^ are the measured photon energies, E[ 2 the re- 
calculated, m T o „ are the nominal ir°,r] masses, and m 77 
the measured invariant masses. 

The obtained 6-photon invariant mass distributions us- 
ing the re-calculated 7-energies are shown in Fig. [5j (left 
column) for the inclusive reaction (np)rf and in coinci- 
dence with recoil protons {n)prj' and neutrons {p)nrf . In 
all cases a clear peak is visible at the nominal ?/-mass of 
958 MeV. The shape of the invariant mass peaks has been 
generated via a Monte Carlo simulation with the GEANT 
package [UJ and fitted to the data together with a back- 
ground polynomial. The peak-to-background is best for 
(n)pr)', intermediate for {np)rf , and worst for {p)nrf . This 
is as expected from the above discussion: events with seven 
'neutral' hits have a much larger chance for combinatorial 
background (for example from i]ir°, rfw , r\n final states, 
when a photon is falsely assigned as neutron or vice versa) 
than events with six 'neutral' and one 'charged' hit. This 
is also reflected in the background of the missing mass 
spectra, which is much more pronounced when the invari- 
ant mass background is not subtracted (compare center 
and right column in Fig. [S]). 

For the missing mass analysis the recoil nucleons were 
treated as missing particles, no matter if they were de- 
tected or not. The missing mass Am of the reaction was 
calculated for quasi-free production of rf mesons off nu- 
cleons via: 



Am 



N 



P v > I - m N , 



(3) 
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where P 7 , Pjv, P^/ are the four- momenta of the incident 
photon, the incident nucleon (assumed to be at rest), and 
the produced r/-meson; mjv is the nucleon mass. The re- 
sulting distributions peak around zero for quasi-free rf 
production and are somewhat broadened by the momen- 
tum distribution of the bound nucleons, which was ne- 
glected. The distributions are shown in Fig. [5j center col- 
umn. They have been constructed for an invariant mass 
window from 930 - 990 MeV (see Fig.[5j left column). Since 
the background cannot be completely removed by cuts on 
invariant mass and missing mass, in the final step the in- 
variant mass spectra have been analyzed (i.e. fitted by 
line shape and background) for each bin of missing mass. 
The resulting missing mass spectra corresponding to the 
invariant mass peaks are shown in Fig. [5j right hand side. 
Background is much reduced and the shapes of the missing 
mass peaks are quite well reproduced by the Monte Carlo 
simulation. The residual background at positive missing 
masses is mainly due to the t/'tt final state, from events 
where the pion escaped detection. It becomes more impor- 
tant at incident photon energies above 1.6 GeV (see also 
Fig. [5]) and extends into the range of the missing mass 
peak. This explains also the deviation of the simulated 
line shape from the data at positive missing mass. There- 
fore at energies above 1.6 GeV only events with missing 
mass between -200 and MeV were accepted. This reduces 
counting statistics by a factor of two but guarantees neg- 
ligible background contamination. The analysis described 
above was done independently for each bin of incident 
photon energy and 77' polar angle. 
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Fig. 6. Time-of-fiight versus energy for charged particles in 
TAPS under the condition of invariant mass signals (cuts on 
invariant mass peaks, no background subtraction) for the 
channel (left hand side), and the rj' channel (right hand side). 
The prominent band corresponds to protons, the less promi- 
nent band to deuterons. 

Absolute cross sections were derived from the target 
density, the incident photon flux, the decay branching ra- 
tio, the detection efficiency for the rj — > 67 decay, and the 
detection efficiency for recoil nucleons. The detection effi- 
ciency was determined via Monte Carlo simulations using 
the GEANT3 package [37] . which included all features of 
the detector system and all software cuts for particle and 
reaction identification. Events for the final states (np)rj' , 
(p)nr)', and (n)prj' were generated evenly distributed in 
phase-space including the effects of nuclear Fermi motion, 
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AE d [MeV] AE^fMeV] 
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Fig. 7. Upper part: missing energy. Left hand side: deuteron 
missing energy, Right hand side: r\ missing energy after cut 
on deuteron missing energy between ±35 MeV. Points: data, 
histograms: simulation of signal and background shapes. Bot- 
tom part: invariant mass spectra. Cut on deuteron missing 
energy ±35 MeV, cut on 77' missing energy: ±25 MeV (left 
hand side), ±50MeV (right hand side). Points: data, fitted with 
background polynomial plus simulated line shape. 



tracked with the GEANT package, and analyzed in the 
same way as the experimental data. The tracking of recoil 
nucleons was done with the GEANT-CALOR program 
package [58| , which is optimized for hadronic interactions 
from the few MeV to several GeV range, including the 
interactions of low energy neutrons. The efficiency correc- 
tion was then done in the usual way as function of incident 
photon energy and meson cm polar angle. 

Finally, an estimate of the cross section for the co- 
herent reaction — > rj' ' d was extracted in the following 
way. Deuterons in the TAPS detector can be identified 
via a time-of-flight versus energy analysis (time-of-flight 
path 1.18 m) as shown in Fig. [HI where only charged hits 
(identified by the veto detectors) are included. The dis- 
tribution at the right hand side of Fig. [5] is obtained for 
6-photon events with an invariant mass corresponding to 
the 77'-meson. The left hand side of the figure shows for 
comparison the result for the rjir channel, which shows 
a much more pronounced deuteron band since it has a 
larger coherent component. Events can then be selected 
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by a cut on the deuteron band, which gives a quite clean 
data sample for the drj' final state. 

However, due to kinematic reasons (mesons from the 
coherent reaction are mostly emitted to forward angles 
due to the nuclear form factor) most deuterons are emit- 
ted into the solid angle of the Barrel, where they cannot 
be identified via time-of-flight. Therefore a more compli- 
cated analysis based on reaction kinematics was necessary. 
The major steps are summarized in Fig. [7] For deuterons 
in the Crystal Barrel, first the Inner detector was used for 
charged particle identification and it was required that the 
deuteron candidate was co-planar with the //-meson. Fur- 
ther exploiting the two-body kinematics of the final state, 
the cm-energies (photon - deuteron center of momentum 
system) derived from the final state four- vectors were com- 
pared to the respective values calculated from the incident 
photon energy for the deuteron (deuteron missing energy) 
and for the 77' (7/ missing energy) . These spectra are com- 
pared in Fig. [7J upper part to Monte Carlo simulations 
of the coherent and breakup process and of residual back- 
ground from r]ir°ir° phase space contributions not related 
to v( production. The corresponding invariant mass spec- 
tra after more or less stringent cuts on missing energy 
are shown in the bottom part of Fig. [71 The final signal 
strength was then extracted via a fit of a background poly- 
nomial and the simulated invariant mass line shape to the 
data. The systematic uncertainty of this procedure was es- 
timated by a variation of the r/' missing energy cut (±50 
MeV, ±25 MeV, -50 MeV - MeV). The first cut cannot 
completely suppress incoherent background, results from 
the most stringent, asymmetrical cut are unfortunately 
limited by the statistical quality of the data. 

3.2 Systematic uncertainties 

Due to the small cross section of //-production, the re- 
quirement of coincident detection of the mesons with re- 
coil particles, and the low trigger efficiency (only trigger 
signals from photons in TAPS) counting statistics were 
low. Therefore statistical uncertainties were rather large 
(on the order of 15 - 35 % for the (p)nr]' final state), 
which makes it difficult to investigate systematic effects 
hidden in statistical fluctuations. Therefore, systematic 
uncertainties were partly extracted from other reaction 
channels, in particular from the 77 —> 67 channel. 

Three different types of systematic uncertainties may 
effect the quasi-free cross sections: overall uncertainties 
which cancel exactly in the comparison of proton and neu- 
tron final state, uncertainties with similar effects, which 
cancel to a large extend in ratios, and uncertainties related 
to specific final states which do not cancel. 

Into the first category fall the systematic uncertainty 
of the incident photon flux, the uncertainty in target thick- 
ness (as well as effects from possible slight misplacements 
of the target), and the uncertainty from the rf decay branch- 
ing ratios. An estimate of the flux uncertainty was ob- 
tained by a comparison of 77-photoproduction [35) cross 
sections obtained by a separate analysis of beam times (A) 
and (C) (see Tab.[JJ. These two beam times used different 



incident electron beam energies, so that the same photon 
beam energies were mapped to different sections on the 
focal plane detector. Furthermore, since (C) was using an 
unpolarized beam and (A) linear polarization, also the en- 
ergy dependence of the flux was different due to the coher- 
ent peak at roughly 1 GeV. Typical deviations between the 
77 cross sections produced from these two beam times are 
at the 5 % level, maximum deviations around 10 %. Since 
the results from both beam times (with approximately 
equal statistical weight) were averaged, we estimate the 
systematic flux uncertainty at 10 %. The overall system- 
atic uncertainty coming from target thickness and posi- 
tioning is on the order of a few per cent. The systematic 
uncertainty of the branching ratios for 77' — > i]ir°ir° —> 67 
is around 7 %. Allowing for some cancellation, we estimate 
a total normalization uncertainty of ~ 15 %. A compari- 
son of the results for quasi-free production off the proton 
to free proton results for both the 77 and the 77' channel 
(see below) did not reveal discrepancies beyond this level. 

The second class of systematic uncertainties is dom- 
inated by the uncertainty of the 77' identification by the 
missing mass and invariant mass analysis and the simu- 
lation of the 77' detection efficiency, which are of course 
related (the better the respective cuts are reflected by the 
simulation the smaller the uncertainty). The simulation of 
the detection efficiency of photons followed by an invariant 
mass analysis for meson identification is very well under 
control for the CB/TAPS setup. This has been tested for 
example via a comparison of the results for the "fp — > prj 
reaction obtained from the analysis of the 77 — > 27 and 
77 — > 67 decay channels 13,15,18 . Agreement with the 
PDG value of the i^^^o / '^^27 ratio is reported in [T5] 
within an uncertainty at the 2 % level. Since errors for 
the 77 — > 27 channel in photon detection or invariant mass 
analysis enter cubed into the 77 — > 3ir° channel, this sets 
stringent limits on the uncertainty. Already a 2 % error in 
photon detection efficiency would result in an 8 % devia- 
tion between the two decay channels. 

A further uncertainty is related to the choice of the 
event generator used for the simulation. Events were gen- 
erated evenly distributed in phase-space, where the ef- 
fects of Fermi motion were modeled using the deuteron 
wave function in momentum space from |52) . Since the 
correction was done as function of incident photon en- 
ergy and cm polar angle of the 77', the angular distribu- 
tion of the 77' mesons itself is not critical. However, de- 
viations could arise if for example final state interaction 
effects modify the correlation between meson polar angle 
and kinetic energy or between meson polar angle and the 
kinematic variables of the recoil nucleon. Possible system- 
atic effects of this kind were investigated with a different 
simulation, where the detection efficiency for ?7'-mesons 
£(2^/, 0^) and the detection efficiency for recoil nucle- 
ons e(T/v,0jv), N=n,p was quasi-factorized and parame- 
terized in dependence of laboratory kinetic energies T, ; < , 
Tjv and polar angles O v i, On of the particles. Typical 
efficiencies arc 10 % for 77' detection (including trigger ef- 
ficiency), 95 % for protons, and 10 - 30 % for neutrons 
(depending on energy and including the identification of 
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the neutron) . These kinematic observables can be directly 
extracted from the measured data (the neutron kinetic 
energy is extracted via the over determined reaction kine- 
matic from the incident photon energy, the measured rj 
four- vector and the measured neutron angles). Therefore, 
an event-by-event efficiency correction with the product 
e{T n i ,O v >) ■ e(Tjv,0jv) becomes possible, which does not 
rely on any model assumptions about the kinematic final 
state variables. This efficiency correction does, however, 
not include the missing mass cut, which depends on in- 
cident photon energy. A correction for this effect was ex- 
tracted from a phase-space simulation. It does not much 
depend on details of the event generator, since it uses only 
the ratio of two different analyses (with and without miss- 
ing mass cut) of the same simulation. Actually as expected 
from Fig.Othe correction factor is close to 2 for a cut from 
-200 MeV to (left half of the peak). The results from 
the two different detection efficiency simulations agreed 
to better than 5 % for all investigated reaction channels, 
and we assume a systematic uncertainty in this range. 

The effects from the background under the invariant 
mass peaks (see Fig. [5]) , which is more important for the 
neutron channel, have not been treated as independent 
systematic uncertainties, they have been included via the 
peak - background separation into the statistical uncer- 
tainties. An additional systematic uncertainty could arise 
from the missing mass analysis. Variations of the accepted 
missing mass range show significant influence on the ex- 
tracted cross section. This, however, does not seem to be 
a problem of the agreement between simulation and data, 
since the shapes of the signals agree well at the left hand 
side of the peaks, but start to disagree at the right hand 
side where background from rj'w final states is expected. 
Therefore, for incident photon energies above 1.6 GcV 
only events with Am <0 have been accepted. However, 
we assign an additional uncertainty rising from 3 % at 
threshold (r/n background starts to contribute above 1.6 
GeV) to 10 % at the maximum energy. Altogether, inde- 
pendent on the reaction channel, we estimate an uncer- 
tainty of 6 % close to threshold up to 12 % at the highest 
incident photon energies, (not including effects of recoil 
nucleon detection). 

The last class of uncertainties are those related to the 
detection of recoil nucleons, which will not cancel in the 
comparison of neutron - proton cross section ratios. The 
detection of the recoil nucleons was included in the simula- 
tions using the GCALOR package [48 , which is optimized 
for this purpose. For the proton, the quality of this simula- 
tion could be cross checked with experimental data for the 
7P — > prj and jp pir°ir° reactions which have been mea- 
sured with the same setup. The proton detection efficiency 
was simply determined as ratio of the number of events 
with detected recoil proton to the total number of events 
from these reactions. The efficiencies have then also been 
simulated and the simulated and measured values agree 
within 8 % for slow protons and 4 % for fast protons. 
Combining all uncertainties except the overall normaliza- 
tion we estimate for the quasi-free proton channel 10 % at 
threshold rising to 15 % at 2.5 GeV. 



For the neutron detection efficiency there are no di- 
rect measurements with the combined TAPS/CB setup 
in Bonn. For the TAPS detector it had been experimen- 
tally determined from the jp —} nir°ir + reaction at the 
MAMI accelerator in Mainz [33]. The results are consis- 
tent with the GCALOR simulation, when the conditions of 
the Mainz setup are used (T„=250 MeV: simulated 18.5 %, 
from data 19.1 %). The neutron detection efficiency of the 
CB was measured at the LEAR ring at CERN [SO]. Re- 
sults from the present GCALOR simulation are in good 
agreement with the LEAR result except for slow neutrons 
(T n < 75 MeV) , where the efficiency is very dependent on 
detector thresholds and the neutron kinetic energy. How- 
ever, in any case it is necessary to determine 'effective' 
neutron efficiencies which take into account the identifi- 
cation of the neutrons out of at least seven neutral hits 
via the invariant mass analysis discussed in Sec. 13.11 This 
could only be done by simulations. The reduction of the 
efficiency under this conditions compared to the situation 
were only neutrons are simulated is substantial, of the or- 
der of 25 % - 35 %. We therefore estimate the absolute 
systematic uncertainty for neutron detection at the 15 % 
level. Altogether a systematic uncertainty of 16 % at re- 
action threshold up to 20 % at highest incident photon 
energies is estimated for the quasi-free neutron channel 
(excluding the overall normalization uncertainty). 

A further systematic uncertainty could arise from the 
misidentification of recoil nucleons. While the loss of events 
is included in the simulated efficiencies, misidentified pro- 
tons might contaminate the neutron sample or vice versa, 
where the first problem is more severe, due to the smaller 
absolute detection efficiency for neutrons. The properties 
of the Inner detector for proton identification have been 
studied in detail in [33] with simulations and data from 
the reaction jp — > ir°p. The main result was that the av- 
erage efficiency for proton detection (somewhat angle de- 
pendent) is 98.9 % (simulation) respectively 98.4 % (data). 
Also determined where the efficiencies of all three layers 
of the detector (from data: 94.8 % (inner layer), 92.9 % 
(middle layer), 88.1 % (outer layer)). Since in this exper- 
iment the condition for neutrons was that no layer had 
responded, only about 0.04 % of protons may be misiden- 
tified as neutrons. The probability that a neutron acti- 
vates a coincidence of two layers (condition for proton) is 
also negligible. The TAPS veto detectors have on average 
an inefficiency for proton detection at the 4 % level, de- 
pending on kinetic energy. However, for TAPS additional 
separation of proton and neutron hits is provided by the 
time-of-flight versus energy analysis. No trace of the typi- 
cal proton band was seen for 'neutral' events and the pos- 
sible contamination of the neutron sample with protons 
could be estimated at the 1 % level. Cross contamination 
of the recoil nucleon samples was therefore negligible. 

The different systematic uncertainties are summarized 
in Tab. [2] It should be noted that the comparison of the 
quasi-free proton data to free proton data, as well as the 
comparison of the two different neutron analyses (see be- 
low) indicate that these estimates are pessimistic. 
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Table 2. Summary of systematic uncertainties for the quasi- 
free reactions. 1 ' photon flux, target thickness, decay branch- 
ing ratios, 2 ' trigger efficiency, rj analysis cuts, rj detection 
efficiency. When two numbers are given first corresponds to 
threshold, second to E~,= 2.5 GeV, and linear interpolation 



source 


•yd — > {n)prj 


7<2 — > {p)nrf 


overall normalization 1 ^ 


15 % 


15 % 


rf detection 2 1 


6 % - 12 % 


6 % - 12 % 


recoil nucleon detection 


8 % - 4 % 


15 % 


total except overall norm. 


10 % - 15 % 


16 % - 20 % 



For the comparison of the quasi-free jp — > 77'p and 
the 777, — > 77'n reactions, systematic uncertainties except 
the ones from the recoil nucleon detection cancel. How- 
ever, these effects can be controlled in an independent 
way. As discussed above, the cross section is constructed 
for 77' mesons in coincidence with recoil protons (cr p ), for 
rj mesons in coincidence with recoil neutrons (c n ), and for 
77' mesons without any condition for recoil nucleons (<r np ). 
Since coherent production processes are very small (see be- 
low), the cross sections must be related by a np = a n + a p . 
Therefore, the neutron cross section can be extracted in 
two independent ways as a n or as a np — a p , one depending 
only on neutron detection efficiency, the other depending 
only on proton detection efficiency. This method has been 
previously tested for 77-photoproduction [35] , where excel- 
lent agreement between the two results was found. Also for 
77' the two results are in good agreement, their weighted 
average (a n ) is given as final result for the neutron cross 
section and the differences between them (shown in Figs. 
I11I13P are an independent estimate of the uncertainties 
introduced by the recoil particle detection. 



4 Reaction Models 

In the absence of any data for polarization observables, 
only a preliminary interpretation of the data in the frame- 
work of reaction models, using model dependent constraints, 
is possible. In 2003 Chiang et al. [30] have developed a 
reggeized model for 77 and 77' production (7/-MAID), which 
they used to analyze the then available proton data. They 
parameterized contributions from nucleon resonances in 
the usual way in terms of Breit-Wigner curves with en- 
ergy dependent widths. Non-resonant Born-terms were ne- 
glected since they were expected to be small at not too 
high photon energies due to the small ^'-nucleon-nucleon 
coupling constant g n 'NN- However, they are included in 
the most recent version of the model used to fit the present 
data. Contributions from i-channel vector meson exchange 
are important and were incorporated via Regge trajecto- 
ries. They found, that already a model including just one 
Sn resonance (W=1960 MeV) together with the Regge 
trajectories could reproduce the available angular distri- 
butions for 77J — > prj'. Some improvement of the fit was 
possible by addition of a P-wave resonance, where a Pn 
or a P13, both around W = 1950 MeV, gave equally good 



results. However, in the meantime, the database for the 
proton has been much improved and, neither the absolute 
magnitude nor the extreme forward - backward asymme- 
try of the early angular distributions [20] have been sup- 
ported by the later experiments [SniTWTg] . so that these 
fits needed to be updated. For this purpose, the model was 
extended by addition of a D13 resonance. It was then fit- 
ted simultaneously to the free proton data from CLAS and 
ELSA, to the present quasi-free proton data, and to the 
quasi-free neutron data. The effects from Fermi-smearing, 
although not important, were taken into account for the 
quasi-free data sets by folding the model results with the 
nucleon momentum distributions. 

In a different approach Nakayama and Haberzettl [32] 
have analyzed the first CLAS-data [3T]. They extended 
their relativistic meson-exchange model |51] for applica- 
tion to the 775 — >■ prj' reaction by introducing contributions 
from spin-3/2 resonances (the earlier version considered 
only spin 1/2-states) and including energy-dependent res- 
onance widths. In addition to the resonance contributions 
nucleonic s- and u-channel diagrams and, more important, 
mesonic /j-channel contributions (/?, u exchange) are con- 
sidered. However, due to the lack of polarization observ- 
ables, they find also different solutions. In Ref. 32 also 
the p- and w-Regge trajectories have been considered in- 
stead of the t-channel rho and omega meson exchanges to 
describe the CLAS data [31]. There, it is found that the 
Regge description yields similar results as the t-channel 
meson exchange. Therefore, in the present work, we con- 
fine ourselves to their model with t-channel meson ex- 
change. The 'minimum' solution with the smallest number 
of resonances in addition to the nucleonic and mesonic cur- 
rents that gives an acceptable fit includes an Sn (1958) and 
Pu(2104) as well as sub-threshold P l3 (1885) and Di 3 (1823) 
states (which can be considered as non-resonant back- 
grounds). In the following, we call this solution (I), whose 
parameter values are summarized in Table I of 32 . So- 
lution (II) includes a further Di3-state at IF=2084 (and 
different parameter values for the other states as sum- 
marized in Table II of |32j). In further exploratory fits 
(Tables III - V in [35]), more resonances were added, but 
did not improve the fit quality significantly. Therefore, 
in the present work we will only discuss solution (I). A 
first analysis of the present quasi-free proton and neutron 
data was done in the following way. The results of the 
fits to the CLAS data for the free proton (solution (I)) 
of |32| ) have been adopted without any parameter change 
and folded with the momentum distribution of the bound 
proton, using the deuteron wave function in momentum 
space [S3] . They are then compared to the quasi-free pro- 
ton results. For the neutron, only the electromagnetic pho- 
ton - resonance couplings of all states have been varied, 
while all other resonance parameters (position, width, de- 
cay branching ratios) have been taken from the proton 
fit. In a second fit, called solution (la), an additional Sn 
resonance was introduced because the neutron data seems 
to show a broad bump at higher incident photon energies, 
which is not apparent for the proton. In the following, we 
refer to solutions (I) and (la) as the NH model. 
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5 Results 

In the following all quasi-free differential cross sections are 
given in the cm (center-of-momentum) system of the inci- 
dent photon and a target nucleon at rest. Apart from the 
immediate threshold region, such cross sections are only 
moderately smeared out by the effects of nuclear Fermi 
motion and can thus be compared almost directly to the 
corresponding results off free nucleons (see [53] for details). 
The angular distributions have been fitted with Legendre 
polynomials 



where the Ai are expansion coefficients. The phase-space 
factor q*,/k* is also evaluated for the above cm system. 
The total cross sections have been extracted from the lead- 
ing Legendre coefficient Aq of these fits. For some previous 
free proton results [THEE]) only angular distributions but 
no total cross sections or total cross sections extracted by 
integration of the angular distributions have been given. 
For this reason, and in order to treat all data samples in 
a consistent way, also for these data sets total cross sec- 
tions have been extracted in this work from the fits of the 
angular distributions. 



5.1 Quasi-free proton cross section 

We first compare the results from the quasi-free 'yd — > 
(n)pr]' reaction to free proton data. The angular distribu- 
tions are summarized in Fig. [9] The total quasi-free cross 
section off the proton is compared to free proton results 
in Fig. [5] This figure shows also the inclusive quasi-free 
cross section (cr np ) of single rj production off the deuteron 
(no condition on recoil nucleon) and the fully inclusive 
cross section (a x ) including contributions e.g. from rj'ir fi- 
nal states. At the highest incident photon energies roughly 
50 % of the yield comes from such meson pairs. For the 
analysis of the single rf channel these multiple meson pro- 
duction reactions have been eliminated by the condition 
that no further mesons have been seen in the detector 
and by the kinematic constraints discussed in Sec. I3.1l for 
events where additional mesons have escaped detection. 

For the comparison of the quasi-free and free proton 
cross sections, one could fold the free cross section data 
with the momentum distribution of the bound nucleons. 
However, simulations have shown, that for photon energies 
above 1.6 GeV this effect is small compared to the uncer- 
tainty of the data. Therefore we compare the unfolded 
data. Apart from a few energy bins close to threshold, the 
shapes of the angular distributions of quasi-free and free 
proton data are in quite good agreement (see also Fig. rT3] 
for a comparison of the fitted Legendre coefficients). For 
the total cross section, shown in Fig. [5] the agreement 
between the present quasi-free data (blue squares) and, 
in particular, the recent high precision proton data from 
CLAS [T!|] (magenta stars) is excellent. The agreement 



with [18] (open circles) is within the systematic uncer- 
tainties. Altogether, no important nuclear effects were ob- 
served for quasi-free ^'-production off the bound proton. 
Therefore, we expect that quasi-free ?/-photoproduction 
off the bound neutron is a reasonable approximation of 
the free neutron reaction. 
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Fig. 8. Total cross sections for total inclusive yd — > Xr]' 
inclusive quasi- free 7c! — > npr/' (a np ) and quasi-free proton 
yd — > (n)pr)' (<7 P ) cross sections. The quasi-free proton cross 
section is compared to the free proton results from |18] (open 
circles) and [19] (magenta stars) . Bar histogram at the bottom: 
systematic uncertainty of quasi-free proton data excluding the 
overall normalization uncertainty. 



5.2 Inclusive quasi-free cross section 

The most simple approach to estimate the behavior of the 
neutron cross section is the measurement of the inclusive 
yd — > (np)rf reaction, where only the rf is detected, pro- 
duction of further mesons is excluded by the missing mass 
cut, and no conditions for the detection of recoil nucle- 
ons are applied. Since coherent contributions are small, 
the result is the incoherent sum of quasi-free proton and 
quasi-free neutron cross section. 

The advantage of this approach is the comparably good 
statistical quality of data without detection of coincident 
neutrons. The angular distributions for this inclusive reac- 
tion are summarized in Fig. [TU] (the total cross section a np 
is included in Fig. [8]). The angular distributions are com- 
pared to the recent CLAS-data for the free proton |19j . 
scaled up by a factor of two. 



I. Jaegle et al.: Quasi- free photoproduction of ^'-mesons off the deuteron 



11 




1 -0.5 0.5 -1 -0.5 0.5 -1 -0.5 0.5 -1 -0.5 0.5 



cos(0;,) 



Fig. 9. Comparison of quasi-free r{ production off the bound proton ((blue) squares) to the free proton data: (black) open 
circles 18;, (black) open crosses: [31], (magenta) stars: |19] . The numbers given in the figure indicate the bin centers in incident 
photon energy (note: first two bins below free nucleon production threshold). Note: results from [311119] partly not exactly for 
the same energy bins as present results. Closest bins or average of overlapping bins chosen. All uncertainties only statistical. 
Lines: Solid (black): Legendre fits to data present data, dashed (red): solution (I) NH model, dotted (blue): ?/-MAlD. 
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Fig. 10. Angular distributions for the inclusive quasi-free process 'yd —¥ (np)r]' of single 77' production (black triangles). 
(Magenta) stars: free proton results from |19| scaled up by factor of two. Note: results from [19] partly not exactly for the same 
energy bins as present results. Closest bins or average of overlapping bins chosen. All uncertainties only statistical. Full (black) 
lines: Legendre fit of present data, dashed (red) lines: solution (I) of NH model, (blue) dotted lines t/-MAID model. 
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Fig. 11. Total cross section for inclusive a np , proton a p , and 
neutron (er n ) final state. Bar histograms: (a„+a p — a np ), curves 
at bottom: ± systematic uncertainty of <r n . Proton (neutron) 
data points slightly displaced to the left (right) for better read- 
ability of the figure. Curves: fits with reaction models model. 
NH model: solid: solution (I), dotted: solution (la), dashed: 
t/'-MAID. 



Agreement between the two data sets therefore signals 
regions where proton and neutron cross sections are more 
or less identical. This is almost perfectly the case (some- 
times with the exception of the extreme backward angles) 
for incident photon energies above 1.9 GeV. This is the re- 
gion, where the angular distributions are strongly forward 
peaked, which in the models is mainly attributed to the 
contribution of i-channel processes. At incident photon 
energies between 1.6 - 1.9 GeV, the inclusive cross section 
is significantly smaller than twice the proton values, in- 
dicating a region, with a n < <j p , which could be a first 
hint to different resonance contributions. For the lowest 
energy bin (1475 MeV) effects of nuclear Fermi motion be- 
come important. Close to the threshold, energy conserva- 
tion asymmetrically favors nucleon momenta anti-parallel 
to the incoming photon momentum, which results in an 
enhancement of meson backward angles (see [33] for de- 
tails). The data are compared to fits with the NH (solution 
(I)) and MAID model. Shown is the incoherent sum of the 
model results for proton and neutron, where the neutron 
couplings have been fitted to the quasi-free neutron data 
(see next section). Agreement with the data for the other 
versions of the NH model, which are not shown, is similar 
to solution (I). 



5.3 Quasi-free neutron cross section 

As discussed in Sec. [21 the quasi-free neutron cross section 
can be extracted by coincident detection of the recoil neu- 
trons or as difference of the inclusive and quasi-free pro- 
ton cross section. The two methods give similar results, 
the averages of the angular distributions are summarized 
in Fig. [T21 For the fitted Legendre coefficients also the two 
individual data sets are shown in Fig. [13] as an estimate of 
systematic uncertainties. The total cross section is com- 
pared in Fig.[TT]to the proton data and to the results from 
the reaction models. Differences between the two extrac- 
tion methods are indicated by the bar histogram in this 
figure. 
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Fig. 12. Angular distributions for the quasi-free 772 — > nrf re- 
action. Only statistical uncertainties. Solid (black) lines: Leg- 
endre fit to data. Dashed (red) lines: solution (I) of NH model, 
dotted (blue) lines: ?7'-MAlD model. 
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Fig. 13. Coefficients of the Legendre Polynomials for the fitted angular distributions. Left hand column: inclusive reaction 
scaled down by factor of 2. Center column (proton targets): quasi-free data (blue squares); free proton data: open crosses [31] 
(omitted for ^3)), open circles [IS], (magenta) stars [IS]. Right hand column (present quasi-free neutron data): (blue) upward 
triangles from neutron coincidence, (black) downward triangles from difference of inclusive and proton data, (red) circles from 
averaged data. Symbols slightly displaced to the left (right) for upward (downward) triangles to make the plot better readable. 
In all plots solid lines: solution (I) NH model, dashed lines: 7/-MAID; for neutron: (magenta) dotted lines CLAS proton data. 



As expected from the discussion of the inclusive data, 
proton and neutron angular distributions are similar in 
magnitude and shape for photon energies above 1.9 GeV. 
However, at lower energies they are significantly different. 
The high precision CLAS proton data [TO] show a kind of 
double bump structure in the total cross section with shal- 
low maxima around W= 1975 MeV and 1^=2080 MeV (cf. 



Fig. [5]), which at the very limit of statistical significance is 
even reflected in the present proton and inclusive data (cf. 
Fig. [TTj) . The earlier free proton CLAS data [31. may also 
contain such structures. In fact, the model calculation of 
Ref. [52], which fits the differential cross sections data of 
Ref. [3T], has predicted such shallow bump structures in 
the total cross section at about the same two energies. 



I. Jaegle et al.: Quasi- free photoproduction of ^'-mesons off the deuteron 



15 



Table 3. Resonance parameters of the r\ MAID model [30]. Resonance positions M and total widths Ftot in MeV. Resonance 
couplings defined by = v/ ' Pnti'^j (in units of lO~ 3 GeV~ 1 ^ 2 ), where /3 Nv i denotes the -/V* — > Nrf branching ratio and, A™ , 
the helicity amplitude; N — p,n and J = 1/2,3/2. P is the neutron/proton ratio of the electromagnetic widths: R = r in /r ip . 



Coupling constants 
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round nucleon born terms g„'NN 
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; e (p, id), g pN N =2.0, 
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=-0.43 (for details see 


EH). 
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Pu 
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0.77 




0.15 


Pl3 


1920 


100 
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5.0 


0.97 


1.4 


D 13 


2150. 


230. 
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-1.4 
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0.02 



Table 4. Resonance parameters of version (I) of the NH model [32] • Notation as in Tab. [3] Background parameters: g v >NN =0.43, 
g P NN =3.3, k, p nn =6.1, g^NN =10.0, k^nn =0, g pT1 '-y =1.25, g un >^ =0.44. The spin-3/2 resonances, Pi 3 and D13, are sub- 
threshold resonances and, as such, they may be considered as part of the background contribution. Further details of model 
(I) may be found in Table. I of Ref. [32]. Version (la) includes in addition a second S11 resonance at W / =2180 MeV and P=110 
MeV. 



Resonance 



M 
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Al/2 



13/2 
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X3/2 
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For the neutron, the first bump is more clearly visible, 
while the second one is suppressed. 

The fitted coefficients of the Legendre representation 
of the angular distributions are summarized in Fig. 1131 
For the quasi-free neutron the results for the two different 
extraction methods and for the average are shown. Agree- 
ment of the two data sets for the even coefficients (Aq, A2) 
is mostly within statistical uncertainties. Some discrepan- 
cies outside statistical uncertainties are observed for A\ 
and A3. These odd coefficients are very sensitive to the 
extreme forward and backward angles, where due to small 
detection efficiencies uncertainties in the data are more 
important. For the same reason no values are given for A3 
for the earlier CLAS-data [3JJ . Due to the relatively small 
angular coverage the fits were not sensitive to it. 

The excellent agreement between free and quasi-free 
proton data is demonstrated in the center column of the 
figure. Effects from nuclear Fermi motion are mostly in- 
significant. The largest effect results again for the odd 
coefficients (Ax, A3) since the asymmetric preference for 
nucleon momenta anti-parallel to the photon momentum 
induces a false forward-backward asymmetry in the quasi- 
free data. However, the effect is small. In case of the most 
precise recent CLAS data, folding with Fermi motion (not 
shown in the figure) improves slightly the agreement for 
the odd coefficients. For a better comparison of proton and 
neutron data, the new CLAS proton data are included as 
magenta, dotted lines into the pictures of the neutron col- 



umn. The largest difference occurs for the Aq coefficient, 
while the results are quite similar for the Ax, A2, and A3 
coefficients. Only close to threshold there could be some 
systematic deviation. 

The neutron data have been fitted with the NH model 
solution (I) (the other solutions give very similar results) 
and the f/-MAID model. The results are shown in Figs. [Til 
ll"2l [131 For both models, all resonance parameters except 
the electromagnetic neutron couplings were taken over 
from, respectively were dominated by, the fits to the pro- 
ton data. For the NH model also a modified version with 
an additional Sn resonance was tested, which, however 
did not much improve the agreement with the data. The 
parameters of the two models are summarized in Tabs. 
13141 For both models agreement between fit and data is 
less good as for the proton data. For the total cross sec- 
tion, the MAID-fit docs not well reproduce the threshold 
region. None of the fits reproduces the region above 2 GeV. 
For the coefficients of the angular distributions particular 
A2 disagrees with the data, which like for the proton rises 
above 2 GeV while it is flat and very small for the fits. 

In this most basic version both models include apart 
from the background terms an Sn and a Pu resonance 
around W=2 GeV (respectively W=2.1 GeV). It must be 
emphasized again, that in both models these solutions are 
by far not unique. Nevertheless, there seems to be agree- 
ment in so far, as both models need an Sn resonance 
close to threshold in order to fit the sharp rise of the total 
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cross section and a Pn state to explain the shape of the 
angular distributions via an S-P-interference term, which 
essentially has a linear cos((9*,) behavior. Although also 
a P13 state and/or the interference with background am- 
plitudes can give rise to such a behavior. For both states 
there are candidates in the Particle Data Group Review 
[3J, the Su(2090) and the Pn(2100) both one star states 
with not well defined positions and widths. In case of 
the Pn both model analyses result in similar positions, 
widths, and neutron/proton ratio of the electromagnetic 
couplings, although the absolute contribution of this state 
is stronger for the NH model. The Sn has similar posi- 
tions but different widths. In case of the MAID model 
the proton coupling is stronger and there is a negative 
sign between proton and neutron coupling, while for the 
NH model the relative sign is positive and the neutron 
coupling is stronger. Consequently, not even for these two 
'dominant' resonances agreement is found between the two 
models. Obviously, further observables must be measured 
to arrive at better constraints for the model analyses. 

Here, one should also note, that already in the La- 
grangian parameterizations of the background terms dif- 
ferences occur between the two models in both, the struc- 
ture of some Lagrangians, as well as in the numerical val- 
ues of coupling constants. As an example for the vector 
meson currents, the coupling constant for the ^-channel uj- 
exchange is positive for the NH model, while it is negative 
for the MAID model (see Tabs. I4T3")) . On the other hand, 
the corresponding coupling for the p-meson is positive in 
both models. This leads to a destructive interference be- 
tween this two terms in the MAID model, while the inter- 
ference is constructive for the NH model. A more detailed 
analysis of the t-channel background terms is therefore 
also necessary. This might profit from data at higher inci- 
dent photon energies, where this contribution dominates. 
Another example are the baryonic background currents 
where the NH model uses a pure pseudo-vector coupling 
at the NNrf vertex, while MAID uses pseudo-scalar cou- 
pling, again giving rise to a relative sign between the am- 
plitudes of the two models. 



5.4 The coherent —> rj'd reaction 

Coherent photoproduction is important due to its direct 
connection to the iso-scalar parts of electromagnetic tran- 
sition amplitudes. However, due to the dependence on the 
nuclear form factor it is strongly suppressed for heavier 
mesons. The results for the total cross section obtained 
from the analysis discussed in Sec. |3J are summarized in 
Fig. Q3] The typical size of systematic uncertainties can 
be estimated by a comparison of the results from the three 
analyses using more or less stringent cuts for the missing 
energy analysis, which is needed to remove background 
from quasi-free processes. All results are on the order of 
only a few nb, the values from the two more stringent cuts 
{AEr,, between ±25 MeV or between -50 - MeV) are in 
reasonable agreement, those from the most open cut (±50 
MeV) may still include a small background contribution. 
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Fig. 14. Total cross section for the coherent reaction 'yd — s* 
drj '. Different symbols correspond to analyses with different 
analysis cuts: closed (red) points: ±25 MeV r( missing mass, 
open circles: ±50 MeV, open (blue) squares: -50 - MeV. Solid 
line: model prediction using the 7/-MAID multipoles. 

The results from the re-fitted reggeized MAID model 
30 discussed above have been used to model the coher- 
ent reaction, using the parameters summarized in Tab. [3J 
The formalism is based on a standard impulse approxi- 
mation without multiple scattering contributions, i.e. the 
nuclear transition operator is taken as a sum of free single 
nucleon operators sandwiched between the deuteron wave 
functions. The general expression for the nuclear ampli- 
tude reads 

Tmj Mi x 1 = 2 J ^ i>l f (p - ~q) t%, i> Mi (p - ^kj , 

(5) 

where k and q are momenta of the initial photon and 
a final meson. The indices Mf, Mi, and A 7 are respec- 
tively the z-projection of the final and initial deuteron 
spin and the photon polarization index. For the deuteron 
wave function, V'm(p), the momentum-space representa- 
tion 

lMp)= £ £(Uf-M L LM L |lM)u L (|p|) (6) 

L=0,2 M L 

x Y LMl (P)xm-m l 

is taken, with \m being the triplet spin wave function. 
For the spatial part ul{p) the parameterization |54j de- 
rived from the Bonn NN potential model was used. For rf 
photoproduction details of the deuteron wave function are 
more important than for lighter mesons since due to the 
large mass already at threshold large momentum trans- 
fers are involved. This causes a strong dependence on the 
behavior of the nuclear potential at small distances, in 
particular the D-w&ve component of the wave function 
becomes important. The calculations predict, that for ex- 
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ample at photon energies around 1.5 GeV the D- wave con- 
tributes already roughly 50 % to the total cross section. 

The operator for rf photoproduction on a single nu- 
cleon has the well known spin structure in the Pauli ma- 
trix representation 

t^,=K^+L^-a, (7) 

with the spin-flip component L and the non spin- flip K . 

The result corresponding to the parameter set from 
Tab. [3] is compared to the data in Fig.Q3J Given the sim- 
plicity of the impulse approximation, neglecting all two- 
body mechanisms as well as final state interaction effects, 
the agreement is quite good, demonstrating that the rel- 
ative contribution of iso-scalar components is well repre- 
sented in the model. 

The predictions for the angular dependence are very 
sensitive to the assumed resonance and background con- 
tributions. For example, at low incident photon energies, 
where the Sn resonance makes a large contribution, the 
spin-flip amplitude L is dominant giving rise to a forward 
peaking of the cross section. At higher incident photon 
energies, where other resonances and the t-channel back- 
ground dominate, the spin-independent part K is strong. 
This part is proportional to sin(0*/) and thus vanishes at 
forward angles. Due to the statistical limitations of the 
data, it was not possible to extract angular distributions. 
A measurement with better statistical quality is highly 
desirable. 



6 Summary and Conclusions 

We have reported the first measurement of ry'-mesons off 
the deuteron. Both, the quasi-free and the coherent reac- 
tion carry important information about the iso-spin com- 
position of the elementary reaction on the free nucleon. 

For the proton case, it has been demonstrated with 
this experiment, that the quasi- free cross section off the 
bound nucleon is very similar to the free proton cross sec- 
tion. No significant nuclear effects e.g. from FSI processes 
have been observed. At the given level of statistical pre- 
cision of the data, even effects from the momentum dis- 
tribution of the bound nucleons are almost insignificant. 
At low incident photon energies, they cause a small arti- 
ficial forward-backward asymmetry in the angular distri- 
butions (coefficient A\ in the Legendre series). At some- 
what higher energies, folding the free proton cross section 
with the momentum distribution improves slightly (within 
systematic uncertainties) the agreement of free and quasi- 
free data for the ^3 coefficient of the angular distributions. 
Only at very low incident photon energies, also the magni- 
tude of the cross section is affected, but in that range free 
nucleon results are also not available or not precise. Agree- 
ment of the quasi-free data with the most recent measure- 
ment off the free proton from the CLAS-collaboration [115] 
was found at a level much below the estimated system- 
atic uncertainty of the present experiment. Consequently, 
the deuteron can be regarded as very well suited target to 
study the 771 — \ rj'n reaction. 



The quasi-free reaction off the bound neutron has been 
studied in two different ways with different sources of sys- 
tematic uncertainty. In one approach, the Ty'-mesons were 
detected in coincidence with the participant neutrons. In 
the second approach, the cross section obtained for co- 
incident participant protons was subtracted from the in- 
clusive results without condition for recoil nucleons. Since 
the detection efficiency for protons and neutrons is very 
different, a comparison of the two results gives a good 
estimate of the systematic uncertainty for recoil particle 
detection. We had previously found excellent agreement 
for a similar analysis of ^-photoproduction [35] . which 
used the same data set. Also for the present analysis of 
?/-photoproduction good agreement is found, which in- 
dicates, that systematic effects arc well under control. 
Only for the extreme forward and backward angular range, 
some discrepancies remain, which effect mainly the A\ co- 
efficient of the Legendre series for the angular distribu- 
tions (see Fig. ITBI . Altogether, the quality of the present 
jn — > nrj data is mostly limited by counting statistics, 
not so much by the systematic effects originating from 
the complications of a quasi-free reaction. 

Proton and neutron cross sections behave similar at in- 
cident photon energies above 2 GeV, where contributions 
from /j-channcl exchange are important. At lower inci- 
dent photon energies, in particular between 1.6 - 1.9 GeV, 
where the proton cross section peaks, the behavior is dif- 
ferent, which might point to different resonance contri- 
butions, but could also arise from changing interference 
terms between the resonances or between resonances and 
non-resonant background. 

The data have been compared to two different mod- 
els, both with contributions from a similar set of nucleon 
resonances and background terms, in particular from t- 
channel mesonic currents. As already pointed out in |32j 
differential cross sections alone cannot uniquely determine 
the contributing reaction mechanisms. Consequently, in 
the framework of both models different solutions can be 
found. Future measurements of polarization observables 
have to clarify the situation. 

Finally, also a first estimate of the coherent 7c? — >• drj' 
contribution at the level of at most a few nanobarn < 
5 nb for all investigated photon energies is a reasonable 
estimate) has been extracted. This reaction is not only 
important for the iso-spin separation of the elementary 
production amplitudes, but aims also at the study of the rf 
nucleon interaction via FSI effects. The extracted results 
are consistent with an impulse approximation, indicating, 
that the iso-spin composition of the model amplitudes is 
at least not unreasonable and that FSI contributions are 
not substantial. However, the statistical limitation of the 
data is even more important at this low cross section level, 
so that effects beyond impulse approximation could not be 
studied. 
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